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ABSTRACT 

Resonant Switch (RS) model of twin peak high-frequency quasi-periodic oscillations (HF QPOs) 
assumes switch of twin oscillations at a resonant point where frequencies of the upper and lower 
oscillations Vrj and v L become commensurable and the twin oscillations change from one pair of 
the oscillating modes (corresponding to a specific model of HF QPOs) to some other pair due 
to non-linear resonant phenomena. The RS model is used to determine range of allowed values 
of spin a and mass M of the neutron star located at the atoll source 4U 1636-53 where two 
resonant points are observed at frequency ratios vu : vl = 3 : 2, 5 : 4. We consider the standard 
specific models of the twin oscillations based on the orbital and epicyclic geodetical frequencies. 
The resonant points are determined by the energy switch effect exhibited by the vanishing of 
the amplitude difference of the upper and lower oscillations. The predicted ranges of the neutron 
star parameters are strongly dependent on the twin modes applied in the RS model. We demon- 
strate that for some of the oscillatory modes used in the RS model the predicted parameters of 
the neutron star are unacceptable. Among acceptable RS models the most promising are those 
combining the Relativistic Precession and the Total Precession frequency relations or their mod- 
ifications. 
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1 Introduction 

In the Galactic Low Mass X-Ray Binaries (LMXB) containing neutron (quark) stars 
quasiperiodic oscillations (QPOs) of X -ray brightness had been observed at low-(Hz ) 
and high-(kHz) frequencies (see, e.g., van der Klisl 2000. 200(§; Barret et al. 2005 ). 



Since the high frequencies are close to the orbital frequency of the marginally sta- 
ble circular orbit representing the inner edge of Keplerian disks orbiting neutron stars, 
the strong gravity effects must be relevant in explaining HF QPOs. 

The HF QPOs in neutron star systems are often demonstrated as two simultane- 
ously observed pairs of peaks (twin peaks) in the Fourier spectra corresponding to 
oscillations at the upper and lower frequencies (vu,vl). The twin peaks at the up- 
per and lower frequencies (vu, vl) substantially change over time (in one observa- 
tional sequence). Sometimes only one of the frequencies is observed and evolved 



dBelloni et alJl2007l:lBouteher et al Il2010t) . Most of the twin HF QPOs in the so called 
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atoll sources ( van der Klisll2006h have been detected at lower frequencies 600 - 800 Hz 
vs. upper frequencies 900 - 1200 Hz, demonstrating a clustering of the twin HF QPOs 
frequency ratios around 3 : 2. This clustering indicates some analogy to the black 
hole (BH) case w here twin peaks with fixed pair of frequencies at the ratio 3 : 2 a re 



req 

usually observed dAbramowicz et a l. 2003; Belloni et al.l 12007c iTorok et al. 2008a,b; 



Boutelier et al. 2010b . Amplitudes of twin HF QPOs in the neutron star systems are 



usua lly stronger and their coherence times are higher th an those of BH sources (see, 
e.g.. iRemillard & McClintockll2006c iBarret et al.ll2005l) . It is probable that a 3 : 2 
resonance plays a significant role also in the atoll sources containing neutron stars. 
However, this case is much more complicated, as the frequency ratio, al though con- 



centrated around 3 : 2, falls in a much wider range than in the BH systems (Toroketal 
2005c lTorokM2009[ iBoutelieret ai1l2010c iMontero & Zanottill2012h . 



In fact, a multi-peaked distribution in the frequency ratios has been observed, i.e., 
more than one resonance could be realized if a resonant mechanism is involved in gen- 
erating the neutron star HF QPOs. For some atoll neutron star sources the upper and 
lower HF QPO frequencies can be traced along the whole observed range, but the prob- 
ability to detect both QPOs simultaneously increases when the frequency ratio is close 
to ratio of small natural numbers, namely 3 : 2, 4 : 3, 5 : 4 - this has been observed in 
the case of six atoll sources : 4U 1636-53, 4U 1608-52 4U 614-09, 4U 1728-34, 
4U 1820-30, 4U 1735-44 dTorokll2009c iBoutelier et aTll2010b . The analysis of root- 
mean-squared-amplitude evolution in the group of six atoll sources shows that the up- 
per and lower HF QPO amplitudes equal each other and alter their dominanc e while 



passing rational frequency ratios corresponding to the datapoints clustering (Torok 
2009). Such an "energy switch effect" ca n be well exp lained in a framework of non- 
linear resonant orbital models as shown in iHoraket all (12009). Moreover, the analysis 
of the twin peak HF QPO amplitudes in the atoll sources (4U 1636-53, 4U 1608-52, 
4U 1820-30, and 4U 1735-44) indicates a cut-off at resonant radii corresponding to 
the frequency ratios 5 : 4 and 4 : 3 respectively implying a possibility that the accretion 
disk inner edge is located at the innermost resonant r adius rather than at the innermost 
stable circular geodesic (ISCO) dStuchhk et al .11201 lb . 

The evolution of the lower and upper twin HF QPOs frequencies in the atoll sources 
suggests (a very rough) agreement of the data distribution with so called hot spot 
models, especially with the relativistic precession (RP) model prescribing the evolu- 
tion of the upper frequenc y by vu = vk and the lower frequency by vl = vk - v r 
( Stella & Vietril 1999 , 1998 ). In rough agreement with the data are also other mod- 
els based on the assumption of the oscillatory motion of hot spots, or accretion disk 
oscillations, with oscillatory frequencies given by the geodetical orbital and epicyclic 
motion: for example it holds for the modified RP1 model with vu = vg and vl = 
v K - v r ( Bursa 2005b. the total pr ecession (TP) model with vu = v g (or vu = v K ) and 
v L = Vg - Vr dStuchhk et al.ll2007h, the tidal disruptio n (TD) model with vy = vk + v r 
and Vl = vk dCadez et al . 2008; Kostic et al J 120091) . or the warped disk oscillations 
(WD) model with v v = 2v K - v r and v L = 2 (v K - v r ) dKatoll2004 120081) . In all of 
these models the frequency difference vu - vl decreases with increase of the magni- 
tude of the lower and upper frequencies in acco rd with trends given by the observational 
data ( Belloni et alj|2007 : Boutelier et al.ll2O10l) . This qualitative property of the obser- 
vational data excludes the simple model of epicyclic oscillations with vu = vg and 
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v L = v r dUrbanec et al.ll2010bl) that works quite well in the case of HF QPOs in LMXB 
containing black holes (Toro k et al, 2005), or requires substantial modification of the 
model of epicyclic oscillations ( Abram owicz et al.ll201 ll) . 

The vu / vl frequency relations, given by a variety of the relevant frequency-relation 
models mentioned above, can be fitted to the observational data for the atoll sources 
containing neutron (quark) stars, e.g., d ata determined for the atoll source 4U 1636-53 
dBarret et al.ll2005t ITorok et al.ll2Q08allbh . The parameters of the neutron (quark) star 
spacetime can be then determined due to the fitting of the observational data by the 
frequency-relation models. The rotating neutron st ars are described quite well by 
the Hartle-Thome geometry (IHartle & Thorne 1968) characterized by three parame- 
ters: mass M, internal angular momentum / and quadmpole moment Q. It is conve- 
nient to use dimensionless parameters a = J/M 2 (spin) and q = QM/J 2 (dimensionless 
quadmpole moment). In the special case when q ~ 1, the Hartle-Thome geometry 
reduces to the well known and well studied Kerr geometry that is convenient for cal- 
culations in strong gravitational field regime because of high simplicity of relevant for- 
mulae. Near-maximum-mass neutron (quark) star Hartle-Thorne models constructed 
for any given equation of state imply q ~ 1 and the Kerr geometry is applicable quite 
correctly in such situations instead of the Hartle-Thorne geometry ( Toro k et al. 2010i 
Urbanec et al. in preparation). High neutron star masses can be expected in the LMXB 
as a result of accretion and in such a case the use of the simple Kerr geometry is justi- 
fied. 

Quality of the fitting procedure appears to be poor for the atoll s ource 4U 1636-53 
showing resonant frequency ratios 3 : 2 and 5 : 4 (ITorok et al .1120121) . Similar very bad 
fit of observational data to the frequency -relation models was found by lLin et al.l (1201 ll) 
for the atoll sources 4U 1636-53 and Sco X-l: the mass and spin determined by the fit- 
ting procedure are given with very large error for any of the applie d models used in 



the paper, as also for som e other models of t he frequency relati ons dMiller et al. 1998 



Osherovich & Titarchuk 1999; Zhang 2004; Zhang et al. 2006; Mukhopadhvav 2009, 
Shi & Lill2009CT chakrabarti et alj |2009t IShfcoi l[ iMukherTe e & Bhattacharyyaj|2012j). 



The strong disagreement of the data distribution and their fitting by the frequency- 
relation models based on the assumption of the geodesic character of the oscillatory 
frequencies related to proper combinations of the geodetical orbital and epicyclic mo- 
tion caused attempts to find a correction of a non-geodesic origin reflecting some im- 
portant physical ingredients - as e.g., influence of the magnetic field of the neutron star 
dBakala et alJ20120201ofeovaf et alJ20oilStuchlfk & Kolosl2012h. or of t hickness of 
non-slender oscillating tori dRezzolla et all2003t IStraub & Sramkova|[2009 ) that could 
make the fitting procedure much better, as shown in lTorok et al.l (12012b . 

However, it is useful to consider another possibility to improve the fitting proce- 
dures that is not based on the non-geodesic corrections and modifies the present orbital 
resonance models based on the assumption of the geodesic orbital and epicyclic mo- 
tion, or disk oscillations with frequencies related to those of the geodetical motion. 
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2 Resonant switch model of HF QPOs 
in neutron star systems 

We propose a new model of twin peak HF QPOs assuming that the twin oscillatory 
modes creating sequences of the lower and upper HF QPOs can switch at a resonant 
point. According to such a Resonant Switch (RS) model non-linear resonant phenom- 
ena will cause excitation of a new oscillatory mode (or two new oscillatory modes) 
and vanishing of one of the previously acting modes (or both the previous modes), i.e., 
switching from one pair of the oscillatory modes to other pair of them that will be act- 
ing up to the following resonant point. We assume two resonant points at the disk radii 
r out and r m , with observed frequencies v^f-, v? ut and v™, vj°, being in commensurable 
ratios p out = n° m : m° at and p m = n m : m m . These resonant frequencies are determined 
by the energy switch effect (Torok 2009); observations put the restrictions v™ > v°' 



out 



and p m < p° at . In the region covering the resonant point at r out the twin oscillatory 
modes with the upper (lower) frequency are determined by the function v° ut (x; M, a) 
(v° ut (;t:; M, a)). Near the inner resonant point at r; n different oscillatory modes given by 
the frequency functions v™(x; M, a) and v™(x; M, a) occur. We assume all the frequency 
functions to be determined by combinations of the orbital and epicyclic frequencies of 
the geodesic motion in the Kerr backgrounds. Such a simplification is correct with 
high precision for neutron (quark) stars with large masses, close to m aximum allowed 
for a given equation of state. As demonstrated in iTorok et al.l d2010h and Urbanec et 
al. (in preparation), the quadrupole moment of Hartle-Thorne geometry correspond- 
ing to near-extreme masses is very close to the Kerr limit and the orbital and epicyclic 
frequencies are very close to those given by the exact Kerr geometry. 

The frequency functions have to meet the observationally given resonant frequen- 
cies. In the framework of the simple RS model, when two resonant points and two 
pairs of the frequency functions are assumed, this requirement enables determination 
of the parameters of the Kerr background describing with high precision the exterior 
of the neutron (quark) star. The "shooting" of the frequency functions to the reso- 
nant points, giving the neutron star parameters, can be realized efficiently in two steps. 
Independence of the frequency ratio on the mass parameter M implies the conditions 

v° ut (x;M, fl ) : v™\x; M, a) = p ont , (1) 
\§{x;M,a) : v£(x;M,a) = p in (2) 

giving relations for the spin a in terms of the dimensionless radius x and the resonant 
frequency ratio p. They can be expressed in the form a° ut (x, p° ut ) and a ln (x, p m ), or in 
an inverse form x° ut (a, p° ut ) and x m (a, p m ). At the resonant radii the conditions 

v° ut = v ° ut (x; M, a), v£ = v™(x; M, a) (3) 

are satisfied along the functions M™'(fl) and Mj" (a) which are obtained by using 
the functions a out (x, p out ) and a ln (x,p m ). The parameters of the neutron (quark) star 
are then given by the condition 

M™l(a) = M™ n (a). (4) 
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The condition (0]i determines M and a precisely, if the resonant frequencies are de- 
termined precisely. If an error occurs in determination of the resonant frequencies, as 
naturally expected, our method gives related intervals of acceptable values of mass and 
spin parameter of the neutron (quark) star. 

Predictions of the RS model have to be tested. First, the predicted ranges of 
the mass and spin parameters have to be confronted with limits predicted by theoretical 
models of the neutron (quark) star structure. Second, these ranges have to be tested by 
the observational limits on mass and spin given by different phenomena observed in 
X-rays coming from the source, e.g., by the profiled spectral lines. Of course, crucial 
will be the test of fitting the twin peak HF QPO data around the resonant points with 
the results of the RS model. 

Here we restrict ourselves to a rough confrontation of the neutron star theoretical 
predictions on the mass and spin. The other tests are under investigation and will be 
presented in a forthcoming paper. 



3 Orbital and epicyclic frequencies of geodetical 
motion in the Kerr geometry 

The formulae for the vertical epicyclic frequency vg and the radial epicyclic frequency 
v r take in the Kerr spacetime ( describing black holes or exterior of near-maximum mass 



neutr on stars) the form (e.g.. lAliev & Galtsovl 11981 
1998tlTorok&"s tuchlik 20051 



Kato et al. 1998; Stella & Vietri 



V7, 



■ a e v^, 



(5) 



where the Keplerian orbital frequency vk and the related epicyclic frequencies are given 
by the formulae 



a g = 



a r = 



1 

2^ 



(GM\ 



\ 'G / 



1/2 



1 



1 



4a 



i _ 6 + 8fl 

X ;t 3 / 2 



x^l 2 + a 
3 a 2 



2n\GMj x 3 ' 2 + a 



1 



3 a 2 



(6) 

(7) 
(8) 



Here x - r/(GM/c 2 ) is the dimensionless radius, expressed in terms of the gravitational 
radius. 

The Keplerian frequency vk(x, a) is a monotonically decreasing function of the ra- 
dial coordinate for any value of the Kerr geometry spin. The radial epicyclic frequency 
has a global maximum for any Kerr black hole spacetime (0 < a < 1) and also the ver- 
tical epicyclic frequency is not m onotonic if the spin is sufficiently high (see, e.g., 
Kato et alJll998uPerez et alJI 19971) . For the Kerr black-hole spacetimes, the locations 
fi r (a), fig (a) of maxima of the epicyclic frequencies v r , vg are implicitly given by 
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the conditions (Torok & Stuchlfl Eool 



P s (x, a) 
/3 r {x, a) 
Pe{x, a) 



- a } {x,a), 
2 x i/z + a 



where j € {r, 6} . 



1 

x 2 
a 



2 a 



a 

x^ + ?' 

a 2 



(9) 
(10) 
(11) 



For any black hole spin, the extrema of the radial epicyclic frequency % r (a) are lo- 
cated above the marginally stable orbit. The latitudinal epicyclic frequency has extrema 
%e (a) located above the photon (marginally bound or marginally stab le) circular orbit 
in the black hole spacetimes with spin a > 0.748 (0.852, 0.952) dTorok & Stuchlik 
120051) . In the Keplerian disks with the inner boundary at x m ~ x ms , the Umiting 
value a = 0.952 is relevant. Note that in the Kerr naked singularity spacetimes (with 
a > 1) the behaviour of the orbital and epicycli c frequencies is much more complex 



dTorok & Stuchli k 2005; Stuchlik & Schee 2012) and also the related optical phenom- 
ena (as profiled lines) demonstrate strong differences in comparison to those created 
in the black hole or neutron star external spacetimes having a < 1 (iStuchlik & Schee 
2010l 2012 ). However, these differences are important only at radii r < M, and be- 



come irrelevant at radii related t o the exterior of the super spinning quark stars having 
(slightly) a > 1 dLo & Linll201 lb . 



4 Frequency relations of the oscillatory mode pairs 

We concentrate attention on the combinations of the relativistic precession (RP) and 
total precession (TP) frequency relations, defined by the frequency mix for the lower 
frequency, their modifications for the upper frequency, including the higher, "beat", 
harmonics, and the tidal disruption (TD) and warped disk (WD) models. For each 
of the frequency relations under consideration we present the function determining 
the resonant radii x n m (a) given by the upper to lower frequency ratio vu : Vl = n : m 
characterized by the parameter 



= (") 



(12) 



Generally, the condition x mm > jc ms has to be satisfied; the marginally stable orbit x ms , 
assumed to correspond to the inner edge of the Keplerian disks, is implicitly determined 
by the function 

a = a ms (x) = ^ (4 - V3x - 2) . (13) 

4.1 Relativistic precession model 

The RP model is determined by the frequency relation 

(14) 

Vk — v r m 



6 



The resonance condition reads 



a = a K/(K - r \x,p) = ^ [4 - V3jc(1-Prp)-2] (15) 



where 

/^rp = (-^-) =(1- VP) . (16) 

Modification RP1. The upper frequency corresponds to the vertical epicyclic fre- 
quency, the modified frequency relation is determined by 

VS " (17) 



v K - v r m 

The resonance function <j e/(K-r) (jc, p) is properly chosen solution of the equation 

p 2 a 2 g -2pa {\ +a r ) + {\ - a r f = 0. (18) 

We shall not give the function explicitly here because of its complexity. 

Modification RPB. The upper frequency corresponds to the beat frequency, the mod- 
ified frequency relation is determined by 



v K -v r m 

The resonance condition reads 



(19) 



= a^ K - r \x, p) ee ^ [4 - V3*(Wrpb)-2] (20) 



3 

where 



( n-m \ 2 ( y/p - 1 \ 2 
\n + m) \ VP + 1 / 



4.2 Total precession model 

The TP model is determined by the frequency relation 

Vg _ n_ 
v e -v r m 

The resonance condition reads 

a = a^\x,p) S ^ (2(p TP + 2)- 
3(/? T p + 1) I 



(22) 



where 



- yjil ~ Pi?) [Mpip + 1) - 2(2p T p + 1)]| (23) 
Ptp = (— ) =(1- <P) ■ (24) 
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Modification TP1. The upper frequency corresponds to the Keplerian orbital fre- 
quency, the modified frequency relation is determined by 



vg — v r m 

The resonance function a K '^ e ~ r '(x, p) is solution of the equation 

(a g - a r f - 2p (a e + a,-) + p 2 = . 
In the explicit form the resonance condition reads 

a = a K/ «>- r \x,p) = 



(25) 



(26) 



= yfx + 



1 



2 • 3 5 / 6 



x 



A 2 / 3 + B 

A 1/3 



A 1 / 3 



( H3pxW _ ^ 

Va+a 1 / 3 ^ 



B 



where 



A = 6p 2 x 5 + y]36p 4 x 10 - B 3 , 
B = 3 1/3 px 3 [4 + (p-4)x] . 



(27) 



(28) 
(29) 



Modification TPB. The upper frequency corresponds to the beat frequency, the mod- 
ified frequency relation is determined by 



Vg + v r n 
Vg - v r m 



(30) 



The resonance condition reads 

^+r)/<.e-r) {x>p) = 



a — a> 



^x 



2(/?tpb + 2) • 



3(/?tpb + 1) 

- -y/a " Ptpb) [3x(p TPB + 1) - 2(2p T PB + D]l (31) 



where 



Ptpb 



In — m\ 2 
\n + ml 



Vp-i 



(32) 



4.3 Tidal disruption model 

The TD model dKostic et alj|2009h is based on the idea of an orbiting hot spot distorted 
by the influence of the tidal effects of the central black hole or neutron star. It is is 
determined by the frequency relation 

vk + v r n 
v K m ' 

The resonance condition reads 



(33) 



a = a (K+r)/K (x, P )s&4- ^3x( l - p TD ) - 2] 



(34) 
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where 



In - m \ 2 
\ m I 



(i - vp) 2 



(35) 



4.4 Warped disk model 

The warped disk oscillations (WD) model dKatoll2004 120081 1 Wagonel 1 999l) assumes 
the frequency of the disk oscillations to be given by combinations of the (multiples) 
of the Keplerian and epicyclic frequencies. Usually, the inertial-acoustic and g-mode 
oscillations and their resonances could be relevant. We consider as an example the fre- 
quency relation 

2vk — v r n 
2 (v K - v r ) m ' 

The resonance condition reads 



(36) 



_(2K-r)/(2K-2r 



r Hx,p) 3 ^ [4- V3jc(1-/>wd)-2] 



where 



fWD - 



2(« — m) 



2n — m 



2(1- vp) 



Vp 



(37) 



(38) 



Resonant switch model applied to the source 
4U 1636-53 



The observational X-ray data obtained for both the atoll and Z-sources indicate feasi- 
bility of the RS model at least in some of the observed sources. We shall test the atoll 
source 4U 1636-53, where the twin peak HF QPOs span a wide range of frequencies 
crossing the frequency ratios 3 : 2, 4 : 3 and finishing at 5 : 4, just near (or at) the in ner 
edge of the accretion disk dBoutelier et"alll2010t lT5r5k1l2009t ^ Istuchlfk et alhoi lh . In 



the case of the 4U 1636-53 source, behaviour of the observed oscillations indicates 
presence of the resonant effects at disk radii where the frequencies are in the ratios 3 : 2 
and 5 : 4 because of the energy switch effect occurring at these frequency ratios ([Ttirok 
2009); therefore, the source is properly chosen for our test of the RS model. 



5.1 The resonant frequencies 

Using the results of Torokl d2009 ). the resonant frequencies determined by the energy 
switch effect are given in the outer resonant point with frequency ratio vu/vl = 3/2 by 
the frequency intervals 



"uo 



Av out = (970 + 30) Hz , 
Av out = (647 + 20) Hz , 



(39) 
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while at the inner resonant point with frequency ratio vu / vl =5/4 the frequency inter- 
vals are 



v- = 



"LO " 



: Av in = (1180 + 20) Hz, 
Av in = (944+ 16) Hz. 



(40) 



Note that the resonant points determined by using the energy switch effect are in accord 
with observational data points cross ing the lines of constant frequency rat ios 3 : 2 and 
5 : 4 as given in the standard papers (Ba rret et al 1 l2005i Bello ni et al ]|2007l) . 



5.2 Theoretical limits on the mass and spin of neutron 
and quark stars 

It is well known that the neutron star mas s surely cannot exceed the critical value 
of M ~ 3.2 M (IRhoades & Ruffinilll974l) . On the other hand, realistic equations 
of state put limit on the maximal mass of neutron stars around M maX N ~ 2.8 M 
( Postnikov et al. 2010|). In fact, the extre mal maximum M maX N ~ 2.8 M is predicted 
by the field theory dMuller & SerotHl996l) . The limit of M maxN ~ 2.5 M Q is predicted 
by Dirac-Brueckner-Hartree- Fock approach in some special case (Miither et al.ll987 ) 
and the variationa l approaches (Akmal & Pandharipande 1 19971: Akm al et al . 1998) and 
other approaches ( Urbanec et al. 2010al) allow for M m „K ~ 2.25 Mm. On the neutron 
star dimensionless spin the limit of a < q maX N = 0-7 has been recently reported, being 
independent of the equation of state dLo & LirfcOl 1 ). 

For the quark stars the maximal mass is expected somewhat smaller in compari- 
son with the neutron stars because of softer equations of state assumed in modelling 
the quark stars, but m asses around M max Q ~ 2 M are allowed (see, e.g., iGlendenning 



2000; lLo & Linll201 lb However, a substantial difference occurs in the limit on max- 



imal spin, since even slightly superspinning states of quar k stars with a max o > 1, ex- 



ceeding the black hole limit, has been recently reported bv lLo & Lid (120 lib , indepen- 
dently of the details of the equation of state for quark matter. Such a sharp difference 
between the limits on the maximal spin of neutron and quark stars could probably be 
explained by the fact that in quark (strange) stars the strong nuclear force helps in bind- 
ing the stars (or could be the binding force in low mass strange stars), while no such 
effect can be present in the case of neutron stars. 

In the Hartle-Thorne models of rotating neutron stars the spin of the star has to 
be related to its rotational frequency. The rotational frequency of the neutron star at 
the atoll source 4U 1636-53 has been o bserved at f mt _ 580 Hz, or f ml = 290 Hz, if 
we observe doubled radiating structure (Stro hmaver & Markwardtil2002l) . Of course, 
such a rotational frequency is substantially lower in comparison with the mass shedding 
frequency, and the Hartle-Thorne model can be applied quite well, predicting spins 
significantly lower than the maximally allowed spin. 

The theory of neutron star structu re then implies for a w ide variety of realistic 
equations of state the spin in the range dHartle & Thornell968i Urbanec et al. in prepa- 
ration) 

0.1<a<0.4. (41) 
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Fig. 1 : Left panel: possible combinations of mass and spin of neutron star in the atoll source 
4U 1636-53 predicted by the RP and TD model using the scatter of the resonant frequencies at 
each of the two resonant points where RP = v K : ( v K - v r ) = 5 : 4 and TD = (y K + v,.) : v K = 3 : 2. 
The right panel gives the detailed information about the intervals of mass and spin relevant for 
this combination of the two frequency relations assuming the resonant switch model. 

Of course, the upper part of the allowed spin range corresponds to the rotational fre- 
quency / rot = 580 Hz, while the lower part corresponds to / rot = 290 Hz. The related 
restriction on the neutron star (near-extreme) mass reads 

M < 2.5 M . (42) 

Therefore, we shall use these restrictions and make a short comment on the results 
of the RS model in the case of three typical combinations considered in our study. 
We plan for a future work to make a detailed study considering Hartle-Thorne models 
constructed for a large variety of equations of state allowing for sufficiently large mass 
of the 4U 1636-53 neutron (or quark) starQ 

5.3 Ranges of 4U 1636-53 neutron star mass and spin implied 
by the RS model 

We have presented above a variety of the oscillatory mode pairs that deserve consid- 
eration and determined the representative resonance functions a Vu ^ VL (x,p) related to 
the oscillatory pairs. Now we shall test all the possible combinations of the oscillatory 
pairs at the established resonant points and give the related intervals of allowed values 
of the mass and spin of the neutron or quark star. We shall discuss in detail two char- 
acteristic cases of the combinations RP-TD, RP-TP and its modification RP1-TP1. In 
the other cases, we shall give only the resulting allowed intervals of the neutron star pa- 
rameters M and a. Of course, each pair of the frequency relations under consideration 
has to be properly ordered at the resonant points. 

We have made the mass and spin estimates by "shooting" the frequency relations to 
the two resonant points for all the frequency relations presented above. The procedure 

1 Notice that quite recent study of the HF QPOs and line profiles observed in 4U 1636-53 leads the authors 
to claim that comb ining these two effects they are able to predict mass of the neutron star to be ~ 2.4 M Q 
fcannaet al.ll2012h . 
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Fig. 2: RP and TD model fits for the atoll source 4U 1636-53. The gray rectangles indicate 
the scatter of the resonant frequencies at each of the two resonant points. The fits are shown for 
the minimal (a), maximal (b) and central (c) allowed values of neutron star mass and spin implied 
by the resonant switch model. The case (a) corresponds to the left black point on the right 
panel of Fig. □ with v° ul (3 : 2) = 1000 Hz, v° ut (3 : 2) = 667 Hz and v>(5 : 4) = 1160 Hz, 
v^(5 : 4) = 928 Hz, (b) corresponds to the right black point on the right panel of Fig. [TJ with 
y™'(3 : 2) = 940 Hz, v° ut (3 : 2) = 627 Hz and v™(5 : 4) = 1200 Hz, vj"(5 : 4) = 960 Hz, and 
(c) corresponds to the central white point on the right panel of Fig. [TJ with v™'(3 : 2) = 970 Hz, 
y™'(3 : 2) = 647 Hz and vjj(5 :4) = 1 1 80 Hz, v£(5 : 4) = 944 Hz. 



of determining the intervals of mass and spin relevant for the combinations of the fre- 
quency relations is demonstrated in Figs. [TJ— 14] The intervals of the mass and spin 
implied by the procedure of the RS model are summarized in Table[TJ Only combina- 
tions giving a RS solution are presented in Table [TJ where for the frequency relations 
under consideration the relevant resonant points are explicitly given. 

The mass and spin estimates of the RS model related to the data of the 4U 1636-53 
source have to be confronted with restrictions on the neutron star mass and spin implied 
by the theoretical models of neutron stars (or quark stars) and the observed rotational 
frequency of the neutron star at 4U 1636-53 presented above. 
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Fig. 3: Left panel: possible combinations of mass and spin of neutron star in the atoll source 
4U 1636-53 predicted by the RP and TP model using the scatter of the resonant frequencies at 
each of the two resonant points where RP = vk : (vk - v r ) = 5 : 4 and TP = v s : (v s — v r ) = 3 : 2. 
The right panel gives the detailed information about the intervals of mass and spin relevant for 
this combination of the two frequency relations assuming the resonant switch model. 

5.3.1 Combination of the RP and TD models 

For the combination of RP and TD models, the RP model has to be related to the in- 
ner resonant point, while the TD model has to be related to the outer resonant point. 
We have demonstrated that for the frequency scatter at the resonant points, given by 
the Eqs. ( f39l > and ( l40b , the range of allowed values of the mass and spin of the neutron 
star is given by 

0.86 < a < 0.94 , 5.03 < — < 6.27 . (43) 

M Q 

The central point of the frequency ranges implies the central estimates of the neutron 
star parameters to be a = 0.91 and M = 5.67 M . For completeness, we present also 
the frequency dependence of the RP and TD oscillatory modes for the limiting values 
of the spin and mass and for the central values in Fig. [2] Qualitatively, these frequency 
dependencies could be in accord with the observational data. 

Clearly, in the case of the combination of the RP and TD oscillatory modes the RS 
model implies parameters of the neutron star that are totally out of the ranges accepted 
by the recent theory of the structure of the neutron or quark stars both for their spin 
and mass. Therefore, the RP-TD combination can be excluded as a realistic expla- 
nation of the observed data in the atoll source 4U 1636-53 because of high values 
of the predicted neutron star parameters. Nevertheless, it is interesting to notice that 
the RP model fits in the allowed range of mass and spin both resonant points quite well, 
especially for the central point with a = 0.91 and M — 5.7 M , while the TD frequency 
relation is restricted just to the region corresponding to the outer resonant point. 

5.3.2 Combination of the RP and TP models 

For the combination of RP and TP models, the RP model has to be related to the inner 
resonant point, while the TP model has to be related to the outer resonant point. For 
the frequency scatter at the resonant points, given by the Eqs. ( 1391 and d40b . we deduce 
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Fig. 4: RP and TP model fits for the atoll source 4U 1636-53. The gray rectangles indicate 
the scatter of the resonant frequencies at each of the two resonant points. The fits are shown for 
the minimal (a), maximal (b) and central (c) allowed values of neutron star mass and spin implied 
by the resonant switch model. The case (a) corresponds to the left black point on the right 
panel of Fig. [3] with v° ut (3 : 2) = 1000 Hz, v° ut (3 : 2) = 667 Hz and y£(5 : 4) = 1160 Hz, 
v^(5 : 4) = 928 Hz, (b) corresponds to the right black point on the right panel of Fig. [3] with 
v° ut (3 : 2) = 940 Hz, v° ut (3 : 2) = 627 Hz and v£(5 : 4) = 1200Hz, v£(5 : 4) = 960 Hz, and 
(c) corresponds to the central white point on the right panel of Fig.[3]with v™'(3 : 2) = 970 Hz, 
v£ ut (3 : 2) = 647 Hz and v£(5 : 4) = 1 180 Hz, v£(5 : 4) = 944 Hz. 



the range of allowed values of the mass and spin of the neutron star given by 

0.29 < a < 0.70 , 2.25 < — < 3.57 . (44) 

M 

The central point of the frequency ranges implies the central estimates of the neu- 
tron star parameters to be a = 0.54 and M = 2.91 M . In Fig. [4] we present again 
the frequency dependence of the RP and TP oscillatory modes for the limiting values 
of the neutron star spin and mass and for the central values of these parameters. Again, 
these frequency dependencies are in accord with the observational data qualitatively. 
As in the case of TD-RP combination, the RP frequency relation meets both resonant 
points frequency intervals, but only slightly, while the TP relation meets only the outer 
resonant point frequency interval. 



14 



Table 1: Intervals of mass and spin of the neutron star in the atoll source 4U 1636-53 implied 
by the procedure of the RS model. The shaded rows represent those combinations of models that 
give acceptable values of spin and mass. Note that the model WD(3:2) gives identical solution 
as TD(3:2). 



Combination of models spin a mass [M/M ] 



RP(3:2)- 


RP(5:4) 


0.68 - 


0.97 


3.60- 


6.88 


RP(5:4)- 


RP1(3:2) 


0.14- 


0.42 


1.98- 


2.49 


RP(5:4)- 


TP(3:2) 


0.29- 


0.70 


2.25 - 


3.57 


RP(3:2)- 


TP 1(5:4) 


0.27- 


0.74 


2.28- 


4.20 


RP(3:2)- 


TPB(5:4) 


0.18- 


0.65 


2.24- 


3.45 


RP(5:4)- 


TPB(3:2) 


0.90- 


0.94 


5.63 - 


6.27 


RP1(3:2) 


-TP(5:4) 


0.25 - 


0.67 


2.11 - 


2.94 


RP1(3:2) 


-TP(5:4) 


0.988 - 


0.994 


3.64- 


3.85 


RP1(3:2) 


-TP 1(5:4) 


0.10- 


0.34 


1.94- 


2.36 


RP1(3:2) 


-TPB(5:4) 


0.992 - 


0.996 


3.65 - 


3.87 


RP1(5:4) 


-TPB(3:2) 


0.9993 - 


0.9996 


3.47 - 


3.59 


RPB(5:4) 


-TP1(3:2) 


0.18- 


0.67 


2.30- 


4.11 


RPB(5:4) 


-TPB(3:2) 


0.72- 


0.84 


4.52- 


5.56 


TP(3:2)- 


TP 1(5:4) 


0.17- 


0.52 


2.07- 


2.93 


TP(3:2)- 


TPB(5:4) 


0.48- 


0.94 


2.84- 


4.71 


TP1(3:2) 


-TPB(5:4) 


0.08- 


0.38 


2.10- 


2.70 


TD(3 


2)- 


TD(5:4) 


0.00- 


0.69 


2.14- 


4.22 


TD(3 


2)- 


RP(5:4) 


0.86- 


0.94 


5.03- 


6.27 


TD(3 


2)- 


RPB(5:4) 


0.34- 


0.73 


2.75- 


4.43 


TD(3 


2)- 


TP 1(5:4) 


0.70- 


0.74 


3.99- 


4.33 


TD(5 


4)- 


TP1(3:2) 


0.38- 


0.71 


2.89- 


4.48 


TD(5 


4)- 


TPB(3:2) 


0.70- 


0.85 


4.38- 


5.58 


WD(3:2) 


-RP(5:4) 


0.86- 


0.94 


5.03- 


6.27 


WD(3:2) 


-RPB(5:4) 


0.34- 


0.73 


2.75- 


4.43 


WD(5:4) 


-RPB(3:2) 


0.00- 


0.86 


2.56- 


7.02 


WD(3:2) 


-TP 1(5:4) 


0.70- 


0.74 


3.99- 


4.33 


WD(5:4) 


-TP1(3:2) 


0.85- 


0.89 


6.64- 


7.59 


WD(5:4) 


-TPB(3:2) 


0.00- 


0.29 


2.56- 


3.21 


WD(3:2) 


-TD(5:4) 


0.00- 


0.69 


2.14- 


4.22 



The RP-TP combination of the oscillatory modes in the RS model implies parame- 
ters of the neutron star that are quite acceptable in the lower end of the allowed ranges 
of spin and mass. Comparing our results with the allowed ranges of spin (Eq. f4~TT > and 
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mass (Eq.l42b we can see that the RP-TP combination could work for 0.3 < a < 0.4 
and mass 2.25 < M/M < 2.5. 

The physical explanation of the RP-TP model could be relatively very simple. 
A hot spot is oscillating in both vertical and radial directions and its radiation is modi- 
fied by the frequencies vg and vg - v r ; approaching the 3 : 2 resonant point, oscillations 
in the vertical direction are successively damped due to non-linear (e.g., tidal) effects 
and the radial oscillations are enforced. Then the Keplerian and the radial epicyclic 
frequencies become important. 

5.3.3 Combination of the RP1 and TP1 models 

Finally, we discuss one important case of combinations of modifications of the RP and 
TP frequency relations. We do not follow all the details, presenting only the results. 
The TP 1 model has to be related to the inner resonant point, while the RP 1 model has to 
be related to the outer resonant point. For the frequency scatter at the resonant points, 
given by the Eqs. fl39l > and d40j >, we deduce the range of allowed values of the mass and 
spin of the neutron star given by 

0.10<a<0.34, 1.94 < < 2.36. (45) 

M 

The central point of the frequency ranges implies the central estimates of the neutron 
star parameters to be a = 0.23 and M — 2.15 M . 

The RP1-TP1 combination of the oscillatory modes in the RS model implies com- 
plete ranges of parameters of the neutron star that are in a good accord with the theo- 
retical limits on the spin (Eq.HTt and mass (Eq.l42l> of the neutron star in 4U 1636-53. 
This combination could thus be considered as one of the best candidates for explaining 
the observed data of the HF QPOs. 

We can see from Table[TJwhere all the results are presented that some combinations 
of the frequency relations are clearly excluded by the theoretical limits on mass and 
spin of neutron stars, while other give quite reasonable restrictions on the mass and spin 
parameters of the neutron star present in the atoll source 4U 1636-53. We plan to test 
further the acceptable combinations of the frequency relations by fitting the relations 
to the observational data related to the vicinity of the resonant points. Our preliminary 
results indicate that the improvement of the fitting procedure precision could be really 
relevant. We believe that than we are able to fix more precisely the proper combination 
of the frequency relations. 

6 Conclusions 

The RS model has been tested for the atoll source 4U 1636-53 demonstrating possible 
two resonant points in the observed data. For relevant pairs of the oscillatory frequency 
relations the range of allowed values of the mass and dimensionless spin of the neutron 
star are determined giving in some cases acceptable pairs of frequency relations, while 
some other pairs are excluded because of predicting unacceptable values of the spin 
and/or mass of the neutron star at 4U 1636-53. We focused our attention on the test of 
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the frequency relations containing geodetical orbital and epicyclic frequencies or some 
combinations of these frequencies. It should be noted that the cause of the switch of 
the pairs of the oscillatory modes is not necessarily tied to the resonant phenomena 
related to the oscillations governed by the frequencies of the geodetical motion, as 
the switch can be related, e.g., to the influence of the magnetic field of the neutron 
star, or the radiation coming from the surface of the neutron star. Moreover, in some 
sources, e.g., the Circinus X-l, we cannot exclude the possibility of chaotic changes 
of frequency relation pairs in a given fixed interval of observed frequencies. We study 
the resonant phenomena first, leaving other cases to the future studies. 

Generally, the RP-TP combination of the RS model, and its modifications, enable 
acceptable explanation of the observational data for 4U 1636-53 source. This should 
be, however, explicitly tested by fitting procedure applied to the observed twin HF QPO 
sequences related to the resonant points. We plan to make such a test in a future work 
and also to test the RS model in the case of some other atoll (4U 1608-52) or Z (Circi- 
nus X-l) sources, containing a neutron (quark) star, with observational data indicating 
possible existence of two resonant points, and to estimate allowed values of the spin 
and mass of the neutron (quark) stars. 

In the special situations related to accreting neutron stars with near-maximum masses 
the Kerr metric can be well applied in calculating the orbital and epicyclic geodetical 
frequencies, as has been done in the present paper, where the results of the mass and 
spin interval findings are in agreement with the assumption of near-maximum masses 
of the neutron stars. In general situations, when the neutron star mass is not close to its 
maximum value allowed by the equation of state, or when a precise fitting procedure 
of the observational data is needed, the Hartle-Thorne geometry describing rotating 
neutron stars has to be considered and the orbital and epicyclic frequencies reflect- 
ing influence of mass, spin and quadrupole moment of the neutron star has to be used. 
Nevertheless, the role of th e quadrupole mom ent is relevant only very close to the inner 
edge of the accretion disk (To roket al 
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